/Journal of Magnetism and Magnetic Materials (2008) 1-0 



1 



o 
o 



Checkerboard order in the t-J model on the square lattice 

D. PottblancPSEE, q WebeJS p. Mila , M. SigristEl 



b 



Laboratoire de Physique Theorique UMR 5152, C.N.R.S. & Universite de Toulouse, F-31062 Toulouse, France 
Institut Romand de Recherche Numerique en Physique des Materiaux (IRRMA), PPH-Ecublens, CH-1015 Lausanne, Switzerland 
^Institute of Theoretical Physics, Ecole Polytechnique Federate de Lausanne, BSP 720, CH-1015 Lausanne, Switzerland 
d Theoretische Physik, ETH-Hdnggerberg, 8093 Zurich, Switzerland 



1. Introduction 
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With constantly improving resolution of experimental techniques, novel features in the global phase 

c$ . diagram of high-T c cuprate superconductors have emerged. One of the most striking is the observation by 

£h I scanning tunnelling microscopy (STM), in the pseudogap phase of underdoped Bi2Sr2CaCu208 +( 5 |T] and 

Ca2- a; Na a ;Cu02Cl2 single crystals ;2j, of a form of local electronic ordering, with a spatial period close 

f5 ■ to four lattice spacings. These observations raise important theoretical questions about the relevance of 

\ such structures in the framework of strongly correlated models. Here, we analyze the stability of new 

inhomogeneous phases |3|4|5j which can compete in certain conditions with the d-wave superconducting 

1 RVB state [6|. Both mean-field [3] (MF) and numerical Variational Monte Carlo results [5] (VMC) will be 

. summarized here. 
\Q ■ 

2. Renormalized Mean-field Theory 

in ■ 

O ■ We describe the doped antiferromagnet by a t— J model, 

: 



T3 



First, we replace the local constraints of no doubly occupied sites by statistical Gutzwiller weights (see 
below) and use a mean-field decoupling in the particle-hole channel to obtain a self-consistent renormal- 
O '. ized MF hamiltonian [7j , 

O ■ 

H M F =~ t Yl 9ij(4,a C j,^ + h - a ) 



d " - \j E 4(^ C U> + fc-c - | X «| 2 ) , (2) 
(yV 

where Xji = ( c j CT c i,o-) and where we explicitly assume an inhomogeneous solution and hence inhomoge- 
neous Gutzwiller factors, 

J = 4(l-x l )(l-x j ) 

9lj (l-7*)(l-tf)-8x i x j \xi j \ 2 + 16\Xii\ A ' 

t r 4xiXj(l - Xi)(l - Xj) I 



113 L (l - *?)(! - x?) + 8(1 - a^)|xyf + 16| X yf 
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where the local hole densities. The \xij\ terms account for the correlations of the probabilities 

between nearest-neighbor sites || and were not included in the original MF treatment j3] of the inhomo- 
geneous state. Our MF scheme applies to the "normal" phase. In contrast, an additional decoupling in 
the particle-particle channel leads to an RVB superconducting state jE] . 

Interestingly, the parent insulating RVB state can be viewed as a staggered flux state (SFP) |9|lU|llj 
MF solution of the above hamiltonian. We show here that, away from half-filling, the SFP spontaneously 
acquires (small) bond modulations. 

3. Results and Comparison with VMC 

Assuming a 4 x 4 supercell with overall rotational symmetry (see Fig. [TJ, the MF equations have 
been solved for an average doping of 1/8. A typical solution depicted in Fig. [T] exhibits orbital currents 
according to the original SFP pattern. However, small modulations of the bond hoppings (c| a Cj jCr ) and 
bond exchange terms (Si ■ Sj) as well as the on-site densities appear spontaneously. Numerical values of 
the bond (site) parameters for t/J = 3 are given in Tabled (Table EJ). The charge-density- wave (CDW) 
modulation is found to be an order of magnitude smaller than the bond-order-wave (BOW) modulations. 
It was argued recently [4j that this instability is connected to nesting properties of the small ellipsoidal 
Fermi pockets of the SFP. 

The MF results have been tested by an exact (within statistical errors) VMC energy minimization and 
evaluation of the order parameters (on a 16 x 16 lattice with mixed boundary conditions, see Ref. [3] for 
details). Comparative values are shown in Tables^&|21 Interestingly, the VMC simulation also predicts 
a dominant BOW component w.r.t the CDW one. However, an overall reduction of the magnitudes of 
these modulations by an approximate factor of 4 compared to MF is seen. 

How such a state can compete with the d-wave RVB superconductor is still unclear. However, our 
calculations |3|5j have suggested that the BOW SFP is less sensitive than the superconducting state to a 
moderate (e.g. nearest-neighbor) Coulomb repulsion. In this competition, the role of the disorder and/or 
the surface might be crucial. 

Tabic 1 

Numerical values of the spin-spin correlations (first 3 lines) and of the real (next 3 lines) and imaginary (last 3 lines) parts of the bond 
hoppings / c !<t c Jct ) ' n BOW stagered flux state for t/J = 3 and x ~ 1/8 (see Table |21 for exact densities). (1) from Ref. (2) with 
corrections in Gutzwiller factors. Two last columns: average and relative BOW modulations (in %). 



Bond # 


1 


2 


3 


4 


5 


6 


aver. 


% (±) 


MF (1) 


-0.346 


-0.172 


-0.313 


-0.236 


-0.164 


-0.125 


-0.220 


50% 


MF (2) 


-0.258 


-0.221 


-0.310 


-0.194 


-0.170 


-0.184 


-0.219 


32% 


VMC 


-0.2149 


-0.2062 


-0.2128 


-0.1871 


-0.1810 


-0.1750 


-0.1963 


10% 


MF (1) 


0.055 


0.057 


0.119 


0.106 


0.092 


0.073 


0.083 


38% 


MF (2) 


0.057 


0.054 


0.061 


0.081 


0.112 


0.051 


0.069 


44% 


VMC 


0.0706 


0.0742 


0.0840 


0.0833 


0.0855 


0.0825 


0.0797 


9.5% 


MF (1) 


0.045 


0.037 


0.037 


0.031 


0.023 


0.030 


0.034 


32% 


MF (2) 


0.038 


0.038 


0.044 


0.031 


0.011 


0.036 


0.033 


50% 


VMC 


0.0400 


0.0392 


0.0414 


0.0349 


0.0313 


0.0330 


0.0367 


12% 
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Fig. 1. Schematic pattern of the bond distribution. Non-equivalent bonds (labelled from 1 to 6) are indicated by different types of lines whose 
widths qualitatively reflect the bond magnitudes. Arrows indicate the directions of the small charge currents. 

Table 2 

Numerical values of the charge densities on the 3 non-equivalent sites (labelled from the center of the 4x4 pattern) for t/J = 3 and doping 
x ~ f/8. Average density and relative amplitudes (in %) are provided in the last column. 



Site # 


f 


2 


3 


n 


% (±) 


MF (f) 


0.925 


0.864 


0.857 


0.8776 


4% 


MF (2) 


0.901 


0.898 


0.806 


0.8759 


5.5% 


VMC 


0.886(3) 


0.874(5) 


0.863(6) 


0.875 


1.3% 



[3] D. Poilblanc, Phys. Rev. B. 72, 060508(R) (2005) 

[4] C. Li, S. Zhou and Z. Wang, Phys. Rev. B. 73, 060501(R) (2006); for related results based on an instability critirium see also, Z. Wang, 

G. Kotliar and X.F. Wang, Phys. Rev. B. 42, 8690 (1990). 

[5] C. Weber, D. Poilblanc, S. Capponi, F. Mila and C. Jaudet, |cond-mat/0601301| 

[6] P.W. Anderson, Science 235, 1196 (1987). 

[7] F.C. Zhang, C. Gros, T.M. Rice and H. Shiba, Supercond. Sci. Technol. 1, 36 (1988). 

[8] T.C. Hsu, Phys. Rev. B 41 11379 (1990); M. Sigrist, T.M. Rice and F.C. Zhang, Phys. Rev. B 49 12058 (1994). 

[9] I. Affleck and J.B. Marston, Phys. Rev. B. 37, R3774 (1988); J.B. Marston and I. Affleck, ibid. 39, 11538 (1989). 



/ Journal of Magnetism and Magnetic Materials (2008) 1-0 



4 



[10] G. Kotliar, Phys. Rev. B 37, 3664 (1988). 

[11] D. Poilblanc and Y. Hasegawa, Phys. Rev. B 41, 6989 (1990); M.U. Ubbcns and P.A. Lee, Phys. Rev. B 46, 8434 (1992); P.A. Lee, 
N. Nagaosa, T. Ng and X. Wen , Phys. Rev. B 57, 6003 (1998); D.A. Ivanov, Phys. Rev. B. 70, 104503 (2004). 



